Weak Lensing Mass Reconstruction: Flexion vs Shear 
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ABSTRACT 

Weak gravitational lensing has proven to be a powerful tool to map directly 
the distribution of dark matter in the Universe. The technique, currently used, 
relies on the accurate measurement of the gravitational shear that corresponds 
to the first-order distortion of the background galaxy images. More recently, a 
new technique has been introduced that relies on the accurate measurement of 
the gravitational flexion that corresponds to the second-order distortion of the 
background galaxy images. This technique should probe structures on smaller 
scales than that of a shear analysis. The goal of this paper is to compare the 
ability of shear and flexion to reconstruct the dark matter distribution by taking 
into account the dispersion in shear and flexion measurements. Our results show 
that the flexion is less sensitive than shear for constructing the convergence maps 
on scales that are physically feasible for mapping, meaning that flexion alone 
should not be used to do convergence map reconstruction, even on small scales. 

Subject headings: Cosmology : Weak Lensing, Shear, Flexion 

1. Introduction 

Weak gravitational lensing is a powerful tool for mapping the distribution of dark matter 
since it measures the matter distribution directly without the need to make assumptions 
about the way that light traces mass. Most approaches focus on shear, which is the first- 
order distortion of the background galaxy images caused by the bending of light from Large 
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Scale Structure (LSS). Several methods have been developed to reconstruct the projected 



mass distribution from the observed shear field (e.g. 


Kaiser et al. 


1993 


Seitz et al. 


1998 


Bridle et al. 


1998 


Marshall et al. 


2002 


Starck et al. 2006 


)■ 



Recently, weak lensing techniques have been extended to include higher-order distortions 
of background galaxies to improve the constraints on the mass distribution on small scales. 
The measurement of the second-order distortion of the background galaxy images by means 



of the galaxy octopole moments was introduced by Goldberg et al. (2002 ). This second-order 



in image distortions corresponding to a third-order effect in gravitational potential (Bacon 



et al. 2006 ) is responsible for the weakly skewed and arc-like appearance of lensed galaxies 
and is expected to probe variations of the gravitational potential field on smaller scales than 
those accessible by shear analysis alone. Despite the fact that measurement of the octopole 
moments is more complex, their intrinsic dispersion due to the random shapes of galaxies 
is expected to be much smaller than the intrinsic ellipticity dispersion. In |Goldberg et al. 



(2005), the method has been further developed by using the shapelet formalism to estimate 



second-order lensing effect. At the same time, a related approach using the galaxy sextupole 



moments has also been explored (Irwin et al. 2003, 2005, 2006). In Okura et al. (2007), the 



authors suggest a new method called HOLICs to measure the second-order lensing effect 
based on the measurement of the octopole and higer-order moments. In |Goldberg et al. 



(2005), the second-order lensing effect was detected for the first time and the term "flexion" 



has been adopted to describe it. The formalism to reconstruct the projected mass distribution 



from the flexion measurements was introduced by Bacon et al. (2006) for the first time 



This paper is structured as follows. In §2, we review the basis of the weak gravitational 
lensing and the flexion formalism. We then proceed with an introduction to the mass inver- 
sion problem from shear and flexion measurements. In §3, a comparison between shear and 
flexion is conducted in order to compare their ability to reconstruct the convergence map 
in the presence of noise. In §4, we have a discussion about published results on mass map 
reconstruction from flexion. In §5, we conclude in the implications of our results on future 
flexion studies. 
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2. Review of weak gravitational lensing formalism 
2.1. Shear formalism 

2.1.1. The second-order of the gravitational lensing potential if} 

Shear 7i(6 l ) with i — 1,2 is measured from the shapes of galaxies at positions 9 in an 
image. The shear field 7«(6 I ) can be written in terms of the lensing potential ip(6) as (see e.g. 



Bartelmann et al. 1999): 



71 = \{dl-d\)i>, 

72 = <9i<9 2 ^, (1) 

where the partial derivatives di are with respect to Q{. The convergence k(8) can also be 
written in terms of the lensing potential as: 

k = \ [p\ + dl) ill. (2) 

The convergence k corresponds to the projected (normalized) mass distribution. 

The left panel of Fig. [T] shows a simulated convergence map derived from ray-tracing 



through N-body cosmological simulations (Teyssier 2002). The cosmological model is taken 



to be a concordance ACDM model with parameters VLm = 0.3, = 0.7, h = 0.7 and 
Og = 0.9. The simulation contains 256 3 particles with a box size of 160/i -1 Mpc. The 
resulting convergence map covers 2x2 degrees with 512 x 512 pixels and assumes a galaxy 
redshift of 1. The typical standard deviation values of k are thus of the order of a few 
percent. 



2.1.2. Shear inversion problem 



The shear mass inversion problem consists of reconstructing the convergence field k(9) 
from the measured shear field "fi(9) by inverting equations Q and ([2]). There are a number 
of approaches for doing this in the literature, and a comparison between the different local 



inversion methods has been carried out by Seitz et al. (1996). 



To simplify the comparison with flexion, we will use the global shear inversion method 



that is presented in Starck et al. (2006) because a similar formalism exists for flexion. For 



this purpose, we take the Fourier transform of the previous equations and obtain 

7< = Pik, i = 1, 2, 



(3) 
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Fig. 1. — The 3 convergence maps have been smoothed with a Gaussian kernel of 15". Left 
: Simulated convergence map k for a ACDM model. The field is 2° x 2° and is obtained from 
N-body simulations with 256 3 particles for a blocksize of \G0h~ 1 Mpc. Middle: Convergence 
map n n reconstructed from noisy shear measurements corresponding to space-based observa- 
tions (crj = 0.3 and n g = 50 gal/arcmin 2 ). Right: Convergence map K n reconstructed from 
noisy flexion measurements corresponding to space-based observations (erf ~ 0.04 arcsec -1 
and n g = 50 gal/arcmin 2 ). 



where the hat symbol denotes Fourier transforms and 

A(k) 



p - k 2 

A, ^ A, i 



2 



2 



A /i \ 2klK2 , . x 

P2(k) = FTP' (4) 

The ideal shear maps 7$ without noise can then be estimated from the convergence map k. 

- 2 ~ 2 

By noting that P\ + P2 = 1, an estimator of the mass distribution k can easily be 
derived by inversion 

kn = Pill + P2l2. (5) 

The deformation induced by weak gravitational lensing on a single galaxy is very weak 
compared to its intrinsic ellipticity. The lensing signal therefore must be extracted from 
the galaxy image ellipticity by assuming the intrinsic ellipticity is randomly oriented in the 
absence of gravitational lensing. The observed shear j i n is then obtained by averaging over 
a finite number of galaxies and, therefore, is noisy. The relationship between the observed 
data 7i, n ,72,n binned in pixels of area A and the true convergence map k are given by: 

7 i,„ = Pi*K + N?, (6) 

where N± and iV^ are white Gaussian noise with zero mean and standard deviation a n ~ 
°1l where N g = n g A is the average number of galaxies in a pixel (n g is the average 

number of galaxies per area unit and A is the pixel area in the same unit). The rms 
shear dispersion per galaxy a] arises both from measurement errors and the intrinsic shape 
dispersion of galaxies. In this analysis, we will assume a] ~ 0.3 as is approximately found 
for ground-based and space-based weak lensing surveys ( Brainerd e al.|1996 ) . Typical values 



for the galaxy surface density for weak lensing are n g ~ 10 gal/arcmin for ground-based 
surveys and n g ~ 50 gal/arcmin 2 for relatively deep space-based surveys. In presence of 
noise, the estimator of the convergence k is: 

k n = Pllln + P2l2n- (7) 

2.2. Flexion formalism 

2.2.1. The third-order of the gravitational lensing potential 



As for shear, the flexion estimation can be calculated with shapelets (Goldberg et al. 
2005 Bacon et al. 2006 Massey et al. 2007) or by directly measuring the higher-order 
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moments of the galaxy image ( Okura et al.||2007 ). Flexion has two components, T and Q. A 
third-order inversion can be performed to recover the convergence field k(9) from the flexion 
field J 7 or Q. It has been shown by Okura et al. (2007) that measurements of the second 
component of flexion Q is more noisy than the first component J 7 . In what follows we will 
only be interested in J 7 . 

The flexion J-{(9) is derived from the second-order shape of galaxies at positions 9 in 
the image. The flexion field J r j(6 l ) can be written in terms of the lensing potential ip(9) as: 

F2 = \{dt + dld 2 )*. (8) 



2.2.2. Flexion inversion problem 

The flexion mass inversion problem consists of reconstructing the convergence field k{9) 
from the measured flexion field JFi{9) by inverting equations (|2]) and (|8]). The Fourier trans- 
form of the relation T\ = d{K gives: 

T\ = —ikik(k), 

f 2 = -ik 2 k(k). (9) 
Then, an estimator of the convergence k can be estimated: 



%k\ - , N ik 



k = ] ^ 2 Fi(k) + w f V2 Mk)- (10) 

In the same way as for shear maps, a measurement error can be associated with the 
flexion maps. The relations between flexion measurements J~ 2 , n and the convergence 
map k are given by: 

J : i t n = d i K + Nf, (11) 

where and N 2 are the noise contributions with a mean equal to zero and a rms equal to 
<y n = erf/ y/Ng- The flexion measurement error erf is between af = 0.01 arcsec -1 (at z = 0) 



and of = 0.1 arcsec" 1 (at z = 1). We choose of ~ 0.04 arcsec -1 as in Bacon et al. (2006). 
In our study, the distribution of the flexion measurements is assumed Gaussian although it 
is not the case in real data. However, this will not affect the conclusions. 

In presence of noise the estimator of the convergence k is: 

k n = j\ 2 An(k) + pjp^,n(fc). (12) 
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3. Comparison 



Flexion should dominate over shear on small scales (Bacon et al. 2006) since flexion 



effects are higher-order deformation of the gravitational potential. Small scale mass distri- 
butions should therefore be covered with higher fidelity with flexion. But what happens 
when measurement errors are added to the data ? 



3.1. Shear noise properties 

The intrinsic ellipticity and the measurement errors on the shear estimation of back- 
ground galaxies result in an additive Gaussian noise on each shear component (see equation 



|6j). The standard dispersion on the shear measurement is oj ~ 0.3 (Brainerd e al. 1996). 
The noise on the convergence map K n is an additive noise iV 7 : 

k n = k + N" 1 (13) 

where : 

j\H = P x Nx + P 2 N 2 \ (14) 

The noise N in k n is still white, Gaussian and uncorrelated. The noise is not amplified 
by the inversion, but k n can be dominated by noise if N is large, which happens in practice. 

To simulate space observations, a realistic white Gaussian noise has been added to 
simulated shear maps. The reconstructed convergence map is dominated by a white gaussian 
noise (a n = 0.181). The middle panel of Fig. [I] shows the reconstructed convergence map 
smoothed by a Gaussian kernel of 15". The smoothing is used to enable the detection of 
some clusters. 



3.2. Flexion noise properties 

The measurement errors on the flexion estimation of background galaxies result in an 



additive Gaussian noise on each flexion component Ti (see equation 11). The dispersion on 
the flexion measurement that comes essentially from the flexion measurement errors is chosen 
to be af ~ 0.04 arcsec -1 . The noise appears on the convergence map K n as an additive noise 
N: 

k n = k + N T , (15) 

where 
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The flexion measurement errors then result in an additive colored Gaussian noise whose power 
is a function of 1/k. The right panel of the Fig. [I] shows a convergence map recovered from 
simulated flexion measurements T^ n . As expected, the convergence map appears contami- 
nated by a colored Gaussian noise whose power is inversely proportional to the frequency k. 
No cluster is detected despite the fact that the map has been smoothed by a Gaussian kernel. 



3.3. Comparison between shear noise and flexion noise 

In this paper, we compare the ability of flexion and shear to reconstruct the dark 
matter distribution. Since flexion dominates on small scales, we calculate here the scale at 
which flexion becomes dominant over shear. To do this, in Fig. [2j we compare the noise 
power spectrum on convergence map obtained from shear measurements (solid black line) 
to the one obtained from flexion measurements (solid red line). The two solid lines have 
been obtained with realistic values of dispersion for space-based observations. The crossing 
of these two curves gives us the scale at which flexion becomes dominant over shear. As 
expected, the shear noise power spectrum is flat and the flexion noise power spectrum is 
inversely proportional to the frequency k. 

10 4 
10 2 
% 10° 

10" 2 
10" 4 

0.001 0.010 0.100 

k [arcsec] -1 

Fig. 2. — Noise power spectrum on convergence map from shear measurements (in black) 
and from flexion measurements (in red) with realistic dispersions (a] = 0.3 and of = 0.04 
arcsec -1 ) and assuming the galaxy density is the same for shear and flexion measurements 
(n g = 50 gal/arcmin 2 ). 




The relations 14 and 16 can be used to derive the analytic shear and flexion noise power 
spectrum: 



a 



7 



2 



\N"i\ z oc — = ,and 



WY oc • (17) 



Nf k 



If the average number of galaxies in a pixel {N g ) is kept the same between shear and flexion 



measurements, the intersection of the two noise power spectra k T is given by k T = If the 
standard values are used for shear and flexion dispersion (crj = 0.3 and erf = 0.04 arcsec -1 , 
k T = 0.1333 arcsec -1 ) this corresponds to a scale of 7.5 arcsec. Thus, the flexion becomes 
interesting for scales smaller than 7.5 arcsec. To have at least a mean of 1 galaxy per pixel 
(with a pixel size of 7.5 arcsec), the galaxy density should be significantly larger than n g ~ 70 
gal/arcmin 2 . Even if one has this galaxy density, all the pixels will not have a galaxy that 
falls inside, and one will have to deal with the problem of missing data. 



Concerning the poor reconstruction of the convergence map from flexion (right panel 
of Fig. [I]), we note that attention has to be paid to the resolution of our simulation. The 
resolution of the simulation that we use for the study is 14 arcsec but many of the scales of 
interest for flexion are below this resolution. However, even with better resolution simula- 
tions these scales are not reachable when mapping with real data - especially if the missing 
data problem is not resolved. 



4. Discussion about published results on convergence reconstruction from 

flexion 

In the literature, several papers have tried to use the flexion to reconstruct the conver- 
gence map. Here is a discussion about the different studies. 



4.1. Non-parametric convergence map reconstruction 



Reconstructions of convergence maps using flexion measurements were first introduced 
by Bacon et al. (2006). In their paper, a (non-parametric) convergence map is reconstructed 
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from simulated flexion measurements. The simulations have a galaxy density of n. 



60 



gal/arcmin 2 and a reported flexion dispersion of <r;f=0.04 arcsec -1 . However, an error in the 
reconstruction code meant that the true flexion dispersion was cr;f=0.007 arcsec -1 , which 
would only be achievable for the highest signal-to-noise galaxie^j] Fig. [3] compares the 
noise power spectrum on the convergence map obtained from realistic shear measurements 
(solid black line) and from realistic flexion measurements (solid red line). The dashed red 



line gives the result from the optimistic flexion measurements used by Bacon et al. (2006) 
(erf =0.007 arcsec -1 and n g = 60 gal/arcmin 2 ). These values are optimistic since to achieve 
this dispersion the flexion of the highest signal-to-noise galaxies should be measured, which 
lead to a galaxy density significantly smaller than n g = 60 gal/arcmin 2 . Doing so increases 

because of the small number of 



17 



the ratio -^r, but decreases the ratio ^ (see equation 
high signal-to-noise galaxies. At the end, the scale k T should remain almost the same. In 



Bacon et al. (2006), the reconstruction fidelity from flexion measurements is therefore too 



optimistic and the result of the reconstruction should be closer to the right panel of Fig. [TJ 



In Okura et al. (2007), the authors also use a (non-parametric) convergence map recon- 
struction from simulated flexion measurements. But the data are simulated with a rather 
optimistic galaxy density n g = 100 gal/arcmin 2 , and a very optimistic value has been chosen 
for the flexion measurement error af = 0.009 arcsec -1 . The dashed red line of Fig. [4] shows 
the noise power spectrum that should be obtained with this optimistic value of flexion dis- 
persion. The intersection with the solid black line gives the scale below which the flexion 
dominates (k T = 50 arcsec). As shown in Okura et al. (2007), the reconstruction of a binned 



convergence map by combining shear and flexion measurements is interesting with this opti- 
mistic flexion dispersion because the flexion is dominating for scales smaller than 50 arcsec. 
But, using real data with a realistic flexion dispersion, the result of the reconstruction of a 
binned convergence map from flexion measurements should be close to what is shown in the 
right panel of Fig. [TJ 



In Okura et aLj (2008), this convergence map reconstruction method is applied to real 
data (ground-based Subaru data). The flexion dispersion is found to be erf = 0.11245 
arcsec -1 . The galaxy density is very small at n g = 7.75 gal/arcmin 2 (only the 791 brightest 
galaxies have been selected) and the field is 9' x 9' sampled with a grid of 256 x 256 pixels. 
Its 791 galaxies to divide into 65536 pixels, which means that only 1% of pixels have a galaxy 
inside. In the paper, there is no mention about the convergence map reconstruction problem 
from incomplete shear maps given that 99 



Vo of data are missing (see Pires et al. (2009) for 
more details about the missing data problem). No detections should be possible in conver- 



1 Private Communication from Bacon, D. 
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Fig. 3. — Noise power spectrum on the convergence map obtained from realistic shear mea- 
surements corresponding to space-based observations (solid black line). The solid red line 
corresponds to the noise power spectrum on convergence map obtained from the realistic 
flexion measurements reported in Bacon et al. (2006) (crT=0.04 arcsec -1 ) and the dashed 



red line corresponds to the noise power spectrum obtained from the very optimistic flex- 
ion measurements that have been really used incorrectly in Bacon et al. (2006) (cr^=0.007 



arcsec 1 ). We assume the galaxy density is the same for shear and flexion measurements 
and we adopt the optimistic galaxy density of Bacon et al. (2006) (n g = 60 gal/arcmin 2 ). 



gence maps obtained from these flexion measurements. But, it is difficult to characterize the 
noise properties of the convergence maps produced by this method and then to access to the 
significance of the detections. 




0.001 



0.010 0.100 
k [arcsec] - 



Fig. 4. — Noise power spectrum on convergence map obtained from realistic shear mea- 
surements (solid black line) and from realistic flexion measurements (solid red line). The 
dashed red line corresponds to the noise power spectrum on convergence map obtained from 
the optimistic flexion dispersion (crf=0.009 arcsec -1 ) of Okura et al, 2007. We assume the 
galaxy density is the same for shear and flexion measurements and we adopt the optimistic 
galaxy density of Okura et al, 2007 (n g = 100 gal/arcmin 2 ). 



In Leonard et al. (2009), an aperture mass is used to reconstruct the convergence map 
from simulated flexion measurements. The flexion dispersion is taken a^ = 0.1 arcsec -1 and 
the galaxy density n g = 35 gal/arcmin 2 . The authors claim to reconstruct substructures, but 
except for the peak that is detected with more than a 3a detection level, no substructure is 
detected with more than a 2a detection level. It should also be noted that a fair comparison 
with the aperture mass for shear is not carried out in this paper. By consequence, no 
conclusions about the utility of flexion measurements to reconstruct substructures can be 
draw from this study. However, the authors are working on this and a paper will be submitted 
later this yearj^] 



2 Privatc communication from Leonard, A. 
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4.2. Parametric convergence map reconstruction 



In Leonard et al. (2007), a parametric convergence reconstruction is performed to re- 
construct A1689 cluster (from HST ACS space-based data), which is one of the biggest and 
most massive known galaxy clusters. The galaxy density is important (n g = 75 gal/arcmin 2 ) 
because of the magnification effect. The measurements were carried out on stacked images, 
which resulted in better shape measurement accuracy (af = 0.029 arcsec -1 ). In this study, 
the galaxy-galaxy flexion signal has been used to show that foreground galaxies are well-fitted 
by a singular isothermal sphere with a characteristic dispersion <j v . Then, for each confirmed 
foreground galaxy, the dispersion a v ^ is estimated from their flexion effect on background 
galaxies. Therefore, the mass reconstruction is modeled as the sum of the fits obtained for 
each foreground galaxy. This method is rather reliable because it depends on the visible 
distribution of the cluster. Then, it offers a way to include the flexion measurements in 
the reconstruction method. However, the measure of the dispersion a v ^ for each foreground 
galaxy remains very noisy and the reconstruction takes no account for the possible presence 
of dark haloes in the cluster. 



Conclusion 



The aim of this paper is to compare the ability of shear and flexion to reconstruct 
convergence maps. A comparison between shear and flexion, taking into account the noise 
contributions, has been carried out. Using noise simulations, we have shown that flexion 
becomes more interesting than shear on scales smaller than the scale containing one galaxy 
(pixel scale). Consequently, the flexion measurements should not be used alone to reconstruct 
a binned convergence map because the flexion is dominating on scales beyond the pixel scale. 
The literature contains several papers that try to use flexion to reconstruct convergence map 
but, their results are not convincing. 

Nonetheless, flexion has already been detected and can still be used to measure the 



statistical properties of substructures in dark matter halos on very small scales (Bacon et al. 



2009). 



Concerning convergence map reconstruction, it is now clear that flexion should not be 
used alone. However, it does help to add the flexion of galaxies in mass reconstruction from 
shear measurements. The question is, how should shear and flexion be combined for optimal 



results ? In Leonard et al. (2007); Shapiro et al. (2010), the authors propose a way of doing 



this but moving beyond this, more working is still needed to find a Bayesian reconstruction 
method for the inclusion of flexion. 
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